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' Abstract: The Nonlinear Schrodinger (NLS) equation is widely used in everywhere of natural science. 
Various nonlinear excitations of the NLS equation have been found by many methods. However, except 



d 
13 



> 



for the soliton-soliton interactions, it is very difficult to find interaction solutions among different types 
of nonlinear excitations. In this Letter, two equivalent very simple methods, the truncated Painleve anal- 
ysis and the generalized tanh function expansion approaches, are developed to find interaction solutions 
between solitons and any other types of NLS waves. Especially, the soliton-cnoidal wave interaction solu- 
tions are explicitly studied by means of Jacobi elliptic functions and the third type of incomplete elliptic 
integrals. The static and moving optical soliton lattices appear naturally in all mediums described the 
NLS systems. Every peak of the cnoidal wave (or soliton lattice) elastically interacts with usual solitons 
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The soliton and/or solitary wave equations connect rich histories of exactly solvable systems con- 
structed in mathematical, statistical and many-body physics, and powerfully demonstrate the unity of 
nonlinear concepts across disciplines and scales from micro-physics and biology to cosmology [T] . Among 
these equations, the nonlinear Schrodinger (NLS) equation is the most ubiquitous 2 . Originally, the NLS 
equation is derived to describe the envelope dynamics of a quasi-monochromatic plane wave propagating 
in a weakly nonlinear dispersive medium when dissipative processes are negligible (see for instance [3]). 
The NLS equation finds an important application in plasma physics, where it describes electron (Lang- 
muir) waves [4]. The NLS equation in nonlinear optics is also well known to describe self- modulation and 
self- focusing of light in a Kerr- type nonlinear medium [5j . The great current interest in the NLS applica- 
tion was initiated by the prediction of the solitons in nonlinear optical fibers [B] and the concept of the 
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soliton laser [7] . Furthermore, the NLS equation is widely used in ferromagnets with easy-axis anisotropy, 
molecular chains, nonideal Bose gas, nuclear matter, solid state medium, gravity waves, optical lattice, 
Bose- Einstein condensations, and so on [8]. 

The multiple soliton solutions of the NLS equation have been given by many authors via different 
methods, say, the Hirota's bilinear method [5], the Darboux transformations (DT) and/or the Backlund 
transformations (BT) |10j . Using the DT and BT, in principle, one can obtain a new solution from a 
known one. However, in practice, one can only find multiple soliton solutions stemming from simple 
constant solutions. It is rather difficult to find new explicit solutions starting from nonconstant nonlinear 
waves such as the cnoidal waves and Painleve waves via DT and BT. In Refs. [TTJ rjj], the mutisoliton 
complexes on a cnoidal wave background have been studied by the DTs for the multi-component NLS 
equations, the sine-Gordon (SG) equation and the Toda lattice. Recently, it is fount that combining the 
symmetry reduction method and the DT or BT related nonlocal svmmeries[13j. one can readily find the 
interaction solutions among solitons and other nonlinear excitations including the cnoidal waves for the 
KdV [H] and KP [15] equations. In this Letter, much more simple and effective methods are proposed 
for the NLS system while they are valid for any integrable systems. 

As an illustration of our new methods, we consider the NLS equation as the special case of the 
Ablowitz-Kaup-Newell-Segur (AKNS) system 

Pt + - ip 2 q = 0, i= V^T, (la) 

qt - Ijibpxx + iq 2 P = 0, b = ±l, (lb) 

with q being a conjugate of p, i.e., q = p*. 

To find significant interaction solutions of the AKNS system among solitons and other types of 
nonlinear excitations, we directly write down the following theorem: 
Theorem. If w is a solution of 
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is a solution of the AKNS system (TT]) with the consistent conditions for the 'phase' u 

u x = 2b— + A, (4a) 
w x 

u t = 3b^l + 4A^ + bw 2 x + hx 2 - h -{w; x}. (4b) 
uu x w x 11 

The consistent condition u x t = u tx of (@| is nothing but @. 

To prove the theorem, we can use two very simple ways, the truncated Painleve analysis and the 

tanh function expansion method. The Painleve analysis is one of the best approaches to study and solve 
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special solutions for nonlinear physical systems. The function (especially the tanh function) expansion 
method is usually used to find traveling wave solutions of the nonlinear partial differential equations. 
Actually, the latter can be considered as the special case of the truncated Painleve expansion approach. 
Using the standard Painleve analysis, it is straightforward to find the following Lemma: 
Lemma. If 6 is a solution of 
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p = Vb 
q = Vb 

is a solution of the AKNS system (TTJ) with the consistent conditions for u 

u x = 2b^- + A, 
3, 



M , = 36^ + 4A^ + ^A 2 

0Z <Px 2 
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Because both the quantities <pt/4> x and {4>; x} are Mobious transformation 
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invariants, the Schwarzian AKNS (and NLS) equation ([S]) is invariant under the Mobious transformation 



Using the Lemma and the following straightening transformation, 
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(9) 



tanh(w 
the theorem is proved. 

We call the transformation (j9)) as the straightening transformation, because it transforms the single 
soliton solution of the Schwarzian AKNS (0 and then the AKNS ([T]) to a straight line solution w — 
k{yx + ui{)t. It is obvious that the solution expression ([3]) is the generalization of the usual tanh function 
expansion method. Actually, the transformation (|9]) converts the usual truncated Painleve expansion 
approach to the most general extension of the special tanh function expansion method. 

Now, the only left thing is to find the interaction solutions of the Eq. ^ with respect to w or 
equivalently the Eq. ([S]) regarding <fi. The theorem reveals that the single soliton (or solitary wave) 
solution of the AKNS system (fTJ) is only a straightened solution w — kx + Lot of Eq. ([2]). This fact hints 
us that to find the interaction solutions between solitons and other nonlinear excitations, we need only 
to find solutions in the form of w = kox + uot + v where v is a function of x and t. Here, we just write 
down the soliton-cnoidal wave interaction solution for {w, u} equations ((2J and (HJ 



w = kox + ujQt + c E^(sn(kx + cut, m), n + 1, m), 
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A ~,w \ , / m 2 ck 3 , 9 3A 2 ,Wn Wn 
u = uo + A + 26-^ )x + b( + k 2 + — + 4A-^ + 3-| 

26c(fcw - k Q u) ^ / . , \ fco(n+l) \ 

+ fcfc0 (fco+M ^ r (fca;+ ^ m) '^w m j > (iob) 

where {uq, m, n, k±, k 2 } are arbitrary constants, other two arbitrary constants related to the initial 

center positions of the soliton and the cnoidal wave have been removed, and 

9 nm 2 „ 9 
c 2 = n , 8 2 = 1, 
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1 + n 

bkd [c 2 k 3 (c 2 -n)- 2ck k 2 (n + nc 2 - n 2 - 2c 2 ) - 3kk 2 (2nc 2 - c 2 - n 2 ) - Anck^] 
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b6k [c 2 k 3 (c 2 -n) + kk 2 (2nc 2 - c 2 - n 2 ) + 2nckf\ 



uj q = ~bk X+ ' =^-. (11) 

2y / ncko(kc + fco)(fcc — nko)(cko + kc 2 — kn) 

In the solution (fT0|) . sn(z,m) is the usual Jacobi elliptic sine function and E^iC, n, to) is the third type 
of incomplete elliptic integral. 

Correspondingly, for the physical quantity, the strength of the AKNS fields {p, q}, I = pq (i.e., 
/ = \p\ 2 for NLS), we obtain 

, [kc-k (S 2 - l)] 2 l2 , , 2bck 2 SCD . - /w„ u\ , 6fc 2 c 2 

/ = &- v ;J tanh 2 (w) - — — tanh(w) + 6A 2 + + - A + 



(1-S 2 ) 2 v ' (l-S 2 ) 2 y ' \ha h) (1-S 2 ) 2 

^^wwi ^ ^c 2 k 3 (c 2 — n) + kk (2nc 2 — c 2 — n 2 ) + 2nck^ ^ ^kc 2 (2n — 1) + n(2cfc — kn) 
kko 2ncko n(S 2 — 1) 

where w is given by (|10[) . S* = yn + lsn(A;x+a;f, m), C = V n + lcn(fca;+u;i, J7i)„ and 13 = dn(fcir+wi, m). 
The soliton-cnoidal wave interaction structure expressed by (|12[) with (| IQf) is abundant. In the following, 
four types of special soliton-cnoidal waves are presented to show a variety of properties of the dark (gray) 
solitons under the background full of bright periodic cnoidal wave. 

Fig. la exhibits the first type of special soliton-cnoidal wave structure of I determined by Eq. ([12")) 
for the NLS equation at t — with the parameter selected as 

{b, ki, k 2 ,S, A, m, n, c, uj , uj} = {1, 1,-1,-1, 0, 0.999, -0.01, 0.00898989, -0.0919418, 0.202157}, (13) 

while Fig. lb displays its time evolution. From Fig. lb, two important features are observed: (i) The 
cnoidal wave looks like a moving optical soliton lattice, (ii) The interaction between soliton and everyone 
of the lattice is elastic except for the phase shifts. 

Fig. 2a displays the second typeof soliton-cnoidal wave structure with the parameter fixed as 

{&, ki,k 2 , (5, A, m, n, c, uj ,uj} = {1, 1, 1,-1, 0, 0.9, -0.2, -0.05, 0.779897, 2.2540497}, (14) 

at the initial time, and its time evolution is explored in Fig. 2b. It is seen from Fig. 2a that for small 
m (here to = 0.9) the cnoidal wave is different from the soliton lattice. However, the interaction between 
the soliton and every peak of the cnoidal periodic wave is also elastic with nonzero phase shifts. 

Fig. 3a, Fig. 3b and Fig. 3c are the time evolutional plots for the third type of soliton-cnoidal wave 
interaction structure with the parameter determined as, 

{b, ki, k 2 , 5, A, to, n, c, ui , ^} = {1, 1, 1, -1, 1, 0.999, -0.875, -2.4720451, -3.0730997, -0.858140}, (15) 



4 




at times t = — 10, and 10, respectively. It is remarkable that if the parameter < n + 1 < lis 
fixed smaller and smaller in Figs. 1-3, then the relative depth of the gray soliton lattice is shallower and 
shallower compared with that of the single non-lattice soliton. 

Fig. 4 shows the fourth type of special soliton-conoidal wave interaction solution with the parameters 
set as 

{b, fa, fa, 6, A, m, n, c, oj , w} = {1, 1, 1, -1, 1, 0.999, -0.00273650, -0.00142257, -0.665981, 0}. (16) 

It is clear that the optical soliton lattice (cnoidal wave) can also be static. In Fig. 4b, the moving soliton 
elastically interacts with every static lattice with the same shift phase. 

It is not surprising that when m — > 1 in Eq. (|10[) . the soliton-cnoidal wave interaction solution 
reduces back to the two-soliton solution, whose interaction behaviors are displayed in Fig. 4 with the 
parameter selections 

{b, fa,fa,S, A, to, n, c, w , u} = {1, 1, 1,-1, 0, 1, -0.95, -4.248259, -5.894616, 0.3634499}. (17) 

To sum up, the NLS system is studied by establishing two quite simple but equivalent methods, the 
truncated Painleve expansion approach and the generalized tanh function expansion method, to find in- 
teractions among different nonlinear excitations. Especially, the soliton-cnoidal wave interaction solutions 
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are explicitly expressed by Jacobi elliptic functions and the third type of incomplete elliptic integral. To 
extend the known method to discover new phenomena in nonlinear physics, the key point is the straight- 
ening transformation (jHJ which straightens the single soliton to a straight line solution and transforms the 
truncated Painleve expansion to an extended tanh function expansion method. The methods are valid 
for all integrable systems and or even nonintegrable models because both the truncated Painleve analysis 
and the tanh expansion method can be used to find exact solutions of the partially solvable nonlinear 
models. The theorem (and the Lemma) can be used to find interaction solutions between solitons and 
any other NLS waves. However, for simplicity, only the soliton-conidal wave interaction solutions are 
discussed in detail. It is found that the soliton-cnoidal interaction solutions possess abundant structures. 
The nonlinear mediums described by NLS system themselves possess optical static or moving soliton 
lattice excitations, the cnoidal wave solutions with the modula very close to 1. Every soliton on the 
soliton lattice elastically interacts with traditional solitons with the only phase shift effect. In general, 
the cnoidal waves with smaller modula are periodic solutions which do not have localized structures. 
However, every peak of the cnoidal waves also elastically interacts with solitons with phase shifts. 

The methods and what obtained in this Letter are valid for all the integrable models. The details 
on the methods for other nonlinear systems, other types of interaction solutions, other methods to solve 
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Fig. 3: Evolution of the third type of special soliton-cnoidal structure of the NLS system with the parameter selections 

JT6l at (a) * = -10, (b) t = and t = 10 respectively. 



interaction solutions among different types of nonlinear excitations and so on will be reported in our 
future research work. 
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